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Abstract.  Apparent  acceleration  of  proper  motion  is  one  of  the  observable  manifestations  of  orbital  motion  in  binary  stars. 
Owing  to  the  increasing  accuracy  of  astrometric  measurements,  it  may  also  be  a  method  to  detect  binarity  of  stars.  This 
paper  presents  some  analytical  expressions  for  the  effects  of  binary  motion  on  proper  motions  when  the  orbital  period  is  at 
least  several  times  the  span  of  observations.  We  estimate  orbit  dimensions  and  distances  at  which  low-mass  companions  and 
planets  may  be  detected  around  main-sequence  stars,  using  preliminary  estimates  of  precision  for  the  AMEX,  GAIA  and 
SIM  space  missions. 
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1.  Introduction 

The  publication  of  the  Hipparcos  catalog  (ESA  1997)  has 
drawn  new  attention  to  the  problem  of  the  determination  of 
proper  motions  of  stars  that  are  components  of  binary  or  mul¬ 
tiple  systems.  With  only  3  years  of  observational  data,  it  was 
recognized  that  unmodeled  binary  motion  could  affect  the 
proper  motions  of  many  stars  in  the  Hipparcos  catalog.  Most 
of  the  concern  has  been  with  binaries  without  known  orbits 
that  may  have  periods  of  several  years  to  several  decades.  Yet 
half  of  solar-type  binary  systems  have  semimajor  axes  of  50 
AU  or  more,  implying,  for  solar-mass  components,  periods  of 
a  century  or  more  (Duquennoy  &  Mayor  1991).  Such  systems 
can  create  problems  when  observational  catalogs  made  years 
or  decades  apart  are  combined  (as  in  Tycho  2)  in  an  effort  to 
provide  improved  proper  motions.  Thus,  the  binary-motion 
problem  occurs  on  many  time  scales  and  may  affect,  to  some 
degree,  a  significant  fraction  of  the  data  in  compiled  catalogs. 

On  the  other  hand,  apparent  accelerations  prove  to  be  a 
powerful  instrument  in  discovering  new  binary  systems.  The 
short-term  Hipparcos  proper  motions  may  be  compared  with 
long-term  ground-based  proper  motions  (e.g.,  fromFK5),  and 
new  binaries  with  a  significant  orbital  motion  can  be  revealed 
even  if  the  companion  is  invisible  or  unresolved.  Wielen  et  al. 
(1999)  developed  a  model  that  estimates  the  statistical  signif¬ 
icance  of  such  differences,  and  pointed  out  that  with  particu¬ 
larly  long  and  accurate  series  of  observations,  it  may  be  pos- 
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sible  to  estimate  the  period  and,  furthermore,  the  mass  of  the 
companion.  Nonlinear  motion  of  photocenters  of  a  few  hun¬ 
dred  suspected  binaries  was  in  fact  estimated  by  Gontcharov 
et  al.  (2001)  from  a  collection  of  ground-based  catalogs  and 
Hipparcos. 

In  this  paper  we  assume  that  no  information  is  available 
for  an  astrometric  analysis  of  proper  motions  other  than  the 
data  obtained  with  a  short-term  space  mission.  We  investi¬ 
gate  stellar  paths  that  are  almost,  but  not  quite,  linear.  This 
case  has  received  less  attention  than  astrometric  orbital  solu¬ 
tions,  although  Eisner  &  Kulkarni  (2001)  presented  an  earlier 
analysis.  We  provide  a  simplified  model  of  astrometric  obser¬ 
vation  analysis  to  quantify  the  problem.  We  use  the  results  to 
estimate  the  range  of  detectable  orbits  for  small  brown  dwarf 
and  Jupiter-mass  companions,  for  three  planned  space  mis¬ 
sions,  GAIA,  SIM  and  the  joint  American-German  project 
AMEX. 

2.  Detectability  of  brown  dwarfs  and  planets 

The  inequality  B3,  derived  in  Appendix  B,  is  used  to  esti¬ 
mate  the  distances  and  orbits  of  very  small  brown  dwarfs 
(O.OIMq)  and  Jupiter-mass  planets  (O.OOIMq)  around  main 
sequence  dwarfs  that  could  be  detected  via  the  acceleration 
method  with  the  AMEX,  GAIA  and  SIM  space  astrometry 
satellites.  The  minimum  detectable  semimajor  axis  with  this 
method  is  defined  by  the  constraint  2 T  <  P,  where  T  is  the 
total  duration  of  a  mission,  and  P  is  the  orbital  period.  Since 
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AMEX 


Fig.  1.  Detectable  orbits  of  a  O.OIA/q  companion  (brown  dwarf). 


AMEX 


Fig.  2.  Detectable  orbits  of  a  O.OOIMq  companion  (a  Jupiter). 

the  mission  length  of  all  three  missions  is  5  years,  the  am -in 
is  the  same  irrespective  of  the  distance  to  the  object,  namely, 
3.4,  4.7  and  6.3  AU  for  a  M2V,  G2V  and  AOV  host,  respec¬ 
tively.  A  single  Jupiter  around  a  Sun-like  star,  therefore,  can 
be  detected  by  any  of  the  satellites  at  appropriate  distances 
(in  the  Solar  system,  ajupiter  =  5.2028  AU). 

The  triangular  shapes  in  Figs.  1  to  5  show  the  domains 
of  detectable  brown  dwarf  and  Jupiter-mass  companions,  us¬ 
ing  acceleration  solutions,  for  the  three  missions  in  question. 
The  detection  domains  of  orbit  solutions  lie  below  these  tri¬ 
angles.  We  assumed  the  mission-average  accuracy  measure 
e  in  Eq.  B3  to  be  twice  the  expected  sky-average  la  accu¬ 
racy  of  parallaxes,  taken  from  the  NASA  SMEX  proposal  for 
AMEX,  (ESA-SCI(2000)4)  for  GAIA,  and  (Unwin  2000)  for 
SIM.  The  upper  limit  amax  for  SIM  is  magnitude-  and  spec¬ 
tral  type-independent,  since  with  this  pointing  mission,  the 
nominal  accuracy  of  parallaxes  («  1  //as)  can  be  achieved 
by  picking  suitable  exposure  times.  As  mentioned  in  Ap¬ 
pendix  A,  the  calculated  upper  detectability  limits  approx¬ 
imately  correspond  to  the  95%  confidence  level;  therefore, 
the  smallest  detected  accelerations  will  have  a  probability  of 
0.95.  Such  moderately  reliable  detections  are  of  interest  for 
the  AMEX  mission  that  would  precede  the  more  accurate 
pointing  SIM  mission  and  select  suitable  targets  for  planet 


GAIA 


Fig.  3.  Detectable  orbits  of  a  0.01  .1/  .  companion  (brown  dwarf). 


GAIA 


Fig.  4.  Detectable  orbits  of  a  O.OOlAfg  companion  (a  Jupiter). 

detection.  The  upper  limit  of  detection  at  a  99%  confidence 
is  fs  1.2  times  lower  in  a  than  the  boundaries  in  Figs.  1-5. 

The  figures  also  show  the  approximate  lower  limits  of  de¬ 
tectable  semimajor  axes  for  orbit  solutions,  although  such  so¬ 
lutions  are  not  considered  in  detail  in  this  paper.  The  thin  dot- 
dashed  lines  in  Figs.  1-5  indicate  the  lower  limit  of  approx¬ 
imate  95%-confidence  detection  domains  for  G2V  primary 
stars  with  orbital  solutions.  It  is  obtained  from  the  simple  re¬ 
quirement  that  the  average  astrometric  excursion  should  be 
larger  than  e/0.87  (cf.  Appendix  2;  0.87  is  an  average  pro¬ 
jection  factor).  Thus,  at  a  given  distance,  and  for  a  given 
primary  and  companion  mass,  the  astrometric  sensitivity  in¬ 
creases  with  semimajor  axis,  provided  at  least  one  full  period 
is  observed.  However,  once  only  a  fraction  of  a  period  can  be 
observed,  the  astrometric  sensitivity  decreases  with  increas¬ 
ing  semimajor  axis,  because  the  observed  motion  becomes 
more  linear.  As  the  distance  increases,  the  width  of  the  regime 
of  detectable  semimajor  axes  narrows,  since  the  observational 
scatter  corresponds  to  an  increasing  linear  scale.  Thus  there 
is  a  reduced  range  of  semimajor  axes  that  are  greater  than 
the  scatter  but  less  than  that  for  which  the  orbital  motion  is 
statistically  indistinguishable  from  a  straight  line. 

Eisner  &  Kulkarni  (2001)  investigated  the  long-period 
regime  from  a  different  perspective,  examining  how  the  pa- 
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SIM 


Fig.  5.  Detectable  orbits  of  a  O.OOlAfg  companion  (a  Jupiter). 

rameters  for  a  simple  orbital  fit  evolve  as  the  period  increases 
beyond  the  span  of  the  observations.  They  considered  several 
cases  of  circular  orbits  with  known  inclinations;  the  analy¬ 
sis  in  our  paper  corresponds  most  closely  to  what  Eisner  & 
Kulkarni  would  call  “type  II  errors”  for  either  face-on  or  in¬ 
clined  orbits.  From  geometric  considerations,  they  obtain  a 
measure  of  the  signal  amplitude  required  for  detection  (their 
eq.  25)  that  is  quadratic  in  P/T  for  P  »  T,  where  P  is 
the  orbital  period  and  T  is  the  time  span  of  the  observations. 
A  little  algebra  shows  that  our  upper  limit  for  the  semima¬ 
jor  axis,  a,  in  eq.  B3,  which  is  proportional  to  T /%/distance, 
is  consistent  with  an  astrometric  signature  that  must  increase 
as  {P/T)'2  for  detection  (all  other  things  being  equal).  Eis¬ 
ner  &  Kulkarni  test  their  analytic  formulas  with  numerical 
simulations.  For  the  face-on  and  inclined  circular  orbits  their 
simulation  results  (their  Figs.  6  and  9)  actually  seem  to  in¬ 
dicate  a  required  amplitude  that  increases  approximately  as 
(P/T)1'9.  It  is  unclear  to  what  extent  this  somewhat  more 
optimistic  result  is  affected  by  the  fact  that  they  apparently 
did  not  solve  for  the  period  (curvature)  of  the  orbits  that  they 
generated. 

3.  Conclusions 

AMEX  will  be  able  to  detect  curved  paths  of  solar-type  stars 
(G2V)  due  to  the  presence  of  a  low-mass  brown  dwarf  out  to 
distances  50  pc.  This  spectral  type  is  optimal  for  a  search 
of  binary  brown  dwarfs  with  AMEX.  The  range  of  detectable 
orbits  at  d  =  10  pc  is  5  to  10  AU.  AMEX  will  be  able,  in 
principle,  detect  the  influence  of  Jupiter-mass  planets  around 
G2V  stars  out  to  ps  4  pc,  and  out  to  m  8  pc  around  M2V  stars, 
but  the  range  of  orbit  sizes  is  fairly  small  at  d  =  4  pc,  3  to  4 
AU.  The  actual  prospect  of  such  detections  are  defined  by  the 
rate  of  jupiters  in  the  immediate  Solar  neighborhood. 

GAIA  will  improve  on  the  detectability  limits  roughly  a 
factor  of  15.  Brown  dwarfs  at  low  orbits  will  be  detectable 
out  to  distances  m  700  pc,  and  jupiters  to  «  100  pc,  for  solar- 
type  stars.  Conversely,  these  objects  may  generate  difficul¬ 
ties  in  the  data  reduction,  requiring  more  complex  models. 
The  rate  of  acceleration  detection  will  grow  at  least  a  factor 


3300,  with  respect  to  AMEX.  Quist  (2001)  conducted  exten¬ 
sive  simulations  of  the  detectability  rates  of  binary  stars  for 
GAIA,  using  a  Galaxy  model  and  provisional  data  reduction 
models.  It  follows  from  his  simulations,  that  more  than  half  of 
all  0.01  A/q.  companions  around  stars  brighter  than  M2V  with 
periods  between  15  and  40  years  will  be  safely  detected  with 
GAIA  (confidence  99.87%).  The  cubic  model,  that  includes 
also  the  second  derivative  of  proper  motion,  can  sometimes 
detect  more  companions  at  smaller  orbital  periods,  than  the 
quadratic  model  considered  in  this  paper,  but  it  is  practically 
useless,  since  at  these  smaller  periods,  the  orbital  model  will 
always  have  a  higher  detection  rate. 

SIM  will  have  moderate  chances  to  detect  Jupiter-mass 
planets  around  late-type  dwarfs  out  to  distances  m  700  pc, 
and  excellent  chances  out  to  ps  200  pc.  But  since  there  are 
perhaps  a  hundred  million  M  dwarfs  within  1  kpc,  SIM  will 
be  able  to  observe  only  a  tiny  fraction  of  them.  A  precursor, 
survey-type  mission  is  therefore  desirable,  that  would  help  to 
select  the  most  promising  targets. 
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Appendix  A:  Simplified  development  of  proper 
motion  estimation 

In  this  section  we  develop  a  simple  model  of  how  analysis  of 
astrometric  observations  for  stellar  proper  motion  determina¬ 
tion  can  become  contaminated  by  an  unmodeled  acceleration. 
The  objective  is  to  provide  approximate  expressions  that  will 
allow  us  to  determine  the  order  of  magnitude  of  the  effect  as 
well  as  its  qualitative  nature. 

The  equation  of  motion  of  a  body  can  be  expressed  as  a 
Taylor  series  in  vector  form  as 

P(i)  =P0+V0f+iZf2  +  ---  (Al) 

where  Po  and  Vo  are  the  body’s  position  and  velocity  at  time 
t  =  0,  and  Z  is  the  acceleration.  For  a  star  where  the  ac¬ 
celeration  is  due  to  the  gravitational  attraction  of  a  compan¬ 
ion  (either  seen  or  unseen),  Z  =  (GM / R2)R,  where  G  is 
the  constant  of  gravitation  and  M  and  R  are  the  mass  and 
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distance  of  the  companion.  The  unit  vector  R  points  toward 
the  companion.  If  we  assume  a  constant  acceleration  —  i.e., 
truncate  eq.  A1  after  the  third  term  —  we  are  limited  to  con¬ 
sidering  a  small  segment  of  the  orbit.  This  is  the  problem  we 
wish  to  investigate:  the  parabolic  approximation  constitutes  a 
“weak  curvature”  case. 

Let  the  star’s  motion,  projected  onto  the  plane  of  the  sky, 
be  p(f).  The  plane  of  the  sky  is  orthogonal  to  the  line  of  sight 
unit  vector  n,  so  we  have 

PW  =  P (t)  ~  (P (t)  •  n)n 

=  (P0  +  Vo t  +  ±Z i2)  -  ((Po  +  V0t  +  \Zi2)  •  n)n 
=  P0  +  V0f  +  \ Zt 2  -  (P0  •  n)n  -  (V0  •  n)th 
-i(Z  •  n)  t2n 

=  [Po  -  (Po  •  n)n]  +  [V0  -  (V0  •  n)n ]t 

+  ±[Z-  (Z-n)n]  t2  (A2) 

That  is,  the  star’s  motion  in  the  plane  of  the  sky  can  be 
represented  as 

P(t)  =  Po  +  v0f  +  t2  (A3) 

where  p0  is  Po  projected  onto  the  sky  =  Po  —  (Po  •  n)n 
Vo  is  Vo  projected  onto  the  sky  =  Vo  —  (Vo  •  n)n 
z  is  Z  projected  onto  the  sky  =  Z  —  (Z  •  n)n 

The  vectors  p (t),  p0,  Vo,  and  z  have  no  component  along  the 
line  of  sight  n  and  are  therefore  2-vectors  in  a  coordinate  sys¬ 
tem  on  the  plane  of  the  sky. 

In  standard  current  practice,  a  star’s  motion  is  expressed 
as  p '(t)  =  Po  +  y't.  Here,  v'  represents  the  star’s  proper  mo¬ 
tion  as  conventionally  defined  —  assumed  constant,  hence 
without  subscript.  (In  this  development  we  are  ignoring  terms 
for  the  curvature  of  the  sky  and  are  assuming  that  aberration 
and  parallax  have  already  been  removed  from  the  data.)  It 
is  tempting  to  think  of  eq.  A3  as  a  simple  extension  of  the 
conventional  expression,  carried  to  higher  order.  But  Vo  in 
eq.  A3  does  not  correspond  to  proper  motion.  For  a  gravita¬ 
tionally  bound  binary,  proper  motion  properly  refers  to  the 
projection  on  the  sky  of  the  (constant)  space  velocity  of  the 
center  of  mass  of  the  system.  In  eq.  A3,  Vo  is  the  instanta¬ 
neous  linear  component  of  the  star’s  apparent  motion,  which 
is  the  sum  of  the  proper  motion  of  system  plus  the  projected 
orbital  velocity  of  the  star  (at  t= 0)  around  the  center  of  mass. 
Determination  of  the  true  proper  motion  of  the  system  would 
require  observations  of  the  star  spanning  nearly  an  orbital  pe¬ 
riod,  or  observations  of  both  the  star  and  its  companion  over 
a  shorter  period  together  with  an  estimate  of  their  mass  ratio. 

But  here  we  wish  to  consider  a  series  of  observations 
where  the  orbital  motion  of  the  star  is  not  obvious.  Consider 
an  observation  of  the  star’s  position,  p(,  taken  at  time  with 
a  measurement  error  e,;  (all  vectors  are  now  in  the  plane  of 
the  sky).  Then  p;  =  p0  A  v0tj  +  i-z t2  +  e,.  But  if  we  have  no 
knowledge  of  the  acceleration,  we  will  model  the  star’s  mo¬ 
tion  in  the  conventional  way  as  p'(i)  =  Pq  +v'£  (see  Fig.  Al). 
The  difference  between  the  observation  and  this  incomplete 
model  of  the  star’s  motion  at  time  I  ,  is  then 

e*  =  Po  +  VqL:  +  \z t2  +  e,:  -  p[,  -  v't,  (A4) 


Compared  to  our  incomplete  model,  the  observation  will  ap¬ 
pear  to  be  in  error  by  an  amount  e' .  This  error  is  the  sum  of  the 
random  error  of  observation  and  the  systematic  error  resulting 
from  the  use  of  the  incorrect  model.  Despite  the  fact  that  it  is 
not  a  purely  stochastic  quantity,  it  provides  a  basis  for  investi¬ 
gating  what  would  happen  if  the  incomplete  model  were  fit  to 
a  set  of  N  two-dimensional  observations  using  least  squares. 
Let  us  take  the  common  situation  where  the  measurement  er¬ 
rors  e,  are  assumed  to  all  have  approximately  the  same  mag¬ 
nitude  (e)  and  the  observations  are  therefore  all  given  unit 
weight.  In  such  a  case,  the  quantity  that  would  be  minimized 
is  the  sum  of  e'2  over  all  observations.  The  method  would 
determine  the  position  \)'(]  (at  time  1=0)  and  proper  motion  v' 
that  minimize  e'2.  The  quantities  p(',  and  v'  are  not  partic¬ 
ularly  interesting  in  themselves,  but  the  differences  between 
these  quantities  and  the  corresponding  ones  from  eq.  A3  for 
this  case  are  key  to  further  analysis.  That  is,  we  are  interested 
in  Ap  =  p0  —  Po  and  Av  =  Vo  —  v'.  These  quantities  measure 
in  some  sense  the  “error”  in  the  position-at-epoch  and  proper 
motion  derived  from  the  linear  fit. 

One  can  perform  a  standard  least-squares  analysis  of  this 
case,  which  amounts  to  the  fit  of  a  straight  line  to  a  series  of 
observations  along  a  parabola;  see  Kaplan  (2002).  We  obtain 

Ap  =  yrj-zT2  and  Av  =  -\zT  (A5) 

Since  the  vectors  Ap,  Av,  and  z  are  collinear,  we  can  dis¬ 
pense  with  the  vector  notation  and  use  use  the  symbols  A p, 
Av,  and  z  for  the  vector  magnitudes.  Equation  A5  then  be¬ 
comes 

A  p=A-zT2  and  Av  = -^zT  (A6) 

The  difference  between  the  actual  position  of  the  star  and 
that  computed  from  the  linear  model  at  time  t  is  in  the  direc¬ 
tion  of  z.  The  acceleration  z  is  always  toward  the  companion, 
and  for  z  to  be  considered  essentially  constant,  the  compan¬ 
ion  must  be  sufficiently  far  away  that  neither  its  direction  nor 
distance  change  significantly  over  the  short  orbital  arc  we  are 
considering.  The  magnitude  of  the  difference  between  the  ac¬ 
tual  motion  and  the  linear  model  is 

8  =  Ap  +  A  vt  +  \zt2  =  pj  zT2  —  \zTt  +  \zt2  (A7) 

This  equation  also  holds  for  the  individual  components  of 
the  motion.  The  function  S  has  an  extremum  (d8/dt= 0)  at 
t=T / 2,  in  the  middle  of  the  span  of  observations,  where 
8  =  —zT2/ 24.  At  the  beginning  and  end  of  the  span  of  ob¬ 
servations  (t= 0  and  f=T),  8  =  zT2  / 12.  The  total  range  of  b' 
over  the  time  interval  of  interest  is  thus  zT2 / 8.  The  locus  of 
actual  motion  of  the  star  on  the  sky  crosses  the  least-squares- 
determined  straight  line  where  the  8  function  has  zeros,  at 
t  =  T/2  ±  sjrjvl  =  0.211  T  and  0.789T.  See  Fig.  Al. 

We  will  refer  to  zT2  / 24,  which  is  the  absolute  value  of 
8  at  t  =  T /2,  as  the  amplitude  of  the  modeling  error  from 
the  linear  approximation,  designated  by  h.  The  total  range 
of  8  over  the  time  span  of  interest  is  3h.  In  the  following 
developments  we  will  use  h  as  the  metric  for  determining 
the  sensitivity  of  the  N  observations  to  the  acceleration.  Intu¬ 
itively,  it  would  seem  that  if  h  is  a  few  times  (e) /s/N,  where 
(e)  is  the  mean  error  of  a  single  observation  (unit  weight), 
then  the  observations  should  be  at  least  marginally  sensitive 
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to  the  acceleration.  More  precisely,  if  we  compute  the  ratio 
of  the  post-fit  sum-of-squares  X2  for  the  linear  and  quadratic 
models,  expressed  as  a  function  of  h,  we  can  use  an  F-test  to 
determine  the  significance  of  the  acceleration  term. 

The  post-fit  sum-of-squares  for  the  case  of  the  linear  fit  to 
accelerated  motion  is 

N 

X2  =N(^z2TA)  +  J2el  (AS) 

i=  1 


Here,  the  e,;  are  the  lengths  of  the  2D  error  vectors,  that  is, 
e2  =  e,  •  e,;.  But  (e)  is  defined  such  that  YliLi  ei  =  A(e)2, 
so  eq.  A8  becomes 

A'  =  N  z2T 4  +  <e)2^  (A9) 


and  the  fractional  increase  in  A"2  (the  “extra  sum  of  squares”) 
due  to  the  modeling  error  from  the  linear  approximation  (and 
assuming  no  other  modeling  errors)  would  be 

2  /  ,  x  2 

(A10) 
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where  2  and  (e)  are  expressed  in  the  same  angular  units,  and 
h  is  the  previously  defined  amplitude  of  the  modeling  error, 
zT2 / 24.  Suppose  we  set  h  =  g(e)/'/N.  where  g  is  a  fac¬ 
tor  to  be  determined;  then  the  fractional  increase  in  X1  is 
g2  /N.  The  ratio  F  that  is  subject  to  the  F-test  is  the  frac¬ 
tional  increase  in  A" 2  (=  fractional  increase  in  y2 ,  since  all 
observations  have  a  an  uncertainty  of  (e))  times  the  number 
of  degrees  of  freedom,  v,  in  the  quadratic  fit  (see,  e.g.,  Bev- 
ington  &  Robinson  1992,  eq.  11.50): 


F  = 


AX2 

xyl 


r 

N1 


(All) 


If  N  is  sufficiently  large  that  v  tx  N,  then  F  22  g2.  For 
an  F-test  probability  of  95%,  F  would  have  to  be  about  4 
for  a  wide  range  of  degrees  of  freedom  (.F=4.17  for  //=30, 
3.92  for  //=120,  3.84  for  v=oo).  Thus  we  require  g  =  2,  i.e., 
h  =  2(e) /\/N.  This  result  is  not  surprising,  since  a  2cr  value 
for  a  model  parameter  determined  from  Gaussian-distributed 
data  has  a  95%  chance  of  being  significantly  different  from 
zero.  If  we  wish  a  more  stringent  test,  we  can  always  set  g=3 
(probability  >  99%);  the  value  of  g  can  be  adjusted  according 
to  the  acceptable  ratio  of  false  positives  /  false  negatives  in  the 
results. 


Appendix  B:  Estimating  the  magnitude  of  the 
effect 

This  section  provides  an  assessment  of  the  magnitude  of 
the  unmodeled  acceleration  effect  that  may  be  detected  with 
a  given  mission-average  precision  of  astrometric  measure¬ 
ments.  The  detectable  orbits  have  apparent  loci  with  a  range 
of  curvatures.  The  development  given  above  treats  the  statis¬ 
tics  of  stellar  loci  at  the  weak  curvature  limit  —  those  repre¬ 
senting  a  very  small  part  (a  few  percent)  of  a  binary  orbital 
period.  The  strong  curvature  limit  might  be  plausibly  defined 
by  the  point  at  which  a  reliable  orbital  solution  becomes  fea¬ 
sible.  That  point  will  be  somewhat  arbitrarily  defined  here  as 


0  0.211T  T/2  0.789T  T 
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Fig.  Al.  Representation  of  a  single  component  of  a  star’s  space  mo¬ 
tion  as  a  function  of  time.  The  true,  accelerated  motion  between 
times  0  and  T  is  indicated  by  the  solid  line  (the  acceleration  compo¬ 
nent  is  zx).  Observational  measurements  are  indicated  that  are  uni¬ 
formly  distributed  in  time  and  normally  distributed  about  the  tiue 
path  of  the  star.  The  dashed  line  is  the  least-squares  straight-line  fit 
to  the  motion,  assuming  equally  weighted  observations.  Various  pa¬ 
rameters  of  the  fit  are  shown. 


being  half  an  orbital  period,  although  preliminary  orbital  so¬ 
lutions  are  often  formed  from  observations  spanning  much 
less  time. 

To  assess  the  practical  effect  of  unmodeled  accelerations, 
a  reasonable  approach  is  to  compare,  for  each  candidate  bi¬ 
nary  system,  the  amplitude  of  the  linear-track  modeling  error, 
h ,  to  some  detection  criterion  e  that  we  are  free  to  choose.  The 
results  of  Appendix  A  show  that  a  reasonable  choice  for  e  is 
2{e)/\/JV,  where  (e)  is  the  mean  error  of  a  single  observation 
and  N  is  the  number  of  observations.  (More  correctly,  the  de¬ 
nominator  should  be  s/v,  where  u  is  the  number  of  degrees  of 
freedom  in  the  quadratic  solution,  but  we  are  assuming  that 
N  is  sufficiently  large  that  v  m  N.) 

The  value  of  the  modeling  error  h  is  defined  for  a  specific 
time  period  of  duration  T,  and  we  will  use  the  expression 
h  =  zT2  / 24  from  Appendix  A.  Although  this  expression 
somewhat  underestimates  the  effect  of  orbital  motion  on  the 
data  for  stars  that  traverse  a  significant  fraction  of  their  orbits 
in  time  T,  it  is  the  appropriate  expression  to  use  for  the  weak 
curvature  limit. 

Our  task,  then,  is  to  determine  over  what  range  of  con¬ 
ditions  h  >  e.  Since  h  is  proportional  to  the  star’s  projected 
acceleration  on  the  sky,  2,  we  must  first  relate  z  to  the  mag¬ 
nitude  of  the  true  acceleration  of  the  star  in  3-space,  |Z|.  The 
true  acceleration  can  be  expressed  in  terms  of  the  physical 
parameters  of  the  binary  system,  and  we  can  then  use  what  is 
known  about  the  distribution  of  these  physical  parameters  to 
estimate  the  frequency  of  significant  acceleration  effects  on 
astrometric  data. 

The  magnitude  of  the  acceleration  projected  onto  the 
plane  of  the  sky  is  related  to  the  magnitude  of  the  true  3-D 


424 


Astron.  Nachr./AN  324,  No.  5  (2003) 


acceleration  vector  by  z  =  (sin#)|Z|/d,  where  #  is  the  an¬ 
gle  between  the  direction  of  Z  and  the  line  of  sight,  and  d  is 
the  distance  to  the  star.  The  parallax,  p,  is  1/d,  so  our  h  >  e 
condition  becomes 

7^-j-  (sin  6)  \Z\T'2  p  >  e  (Bl) 

The  acceleration  magnitude  is  related  to  the  physical  state  of 
the  binary  system  through  |Z|  =  GM/R 2,  where  R  is  the 
instantaneous  distance  between  components,  and  M  is  either 
(a)  the  mass  of  the  companion,  if  the  star’s  motion  is  mea¬ 
sured  in  an  inertial  system,  or  (b)  the  total  mass  of  the  system, 
if  the  motion  is  measured  with  respect  to  the  companion.  If 
we  use  units  of  AU,  years,  and  solar  masses,  the  gravitational 
constant  G  =  Air'2  and  the  true  acceleration  |Z|  is  expressed 
in  AU  year-2.  If  d  is  in  parsecs  and  p  in  arcseconds,  then  2  is 
in  units  of  arcsec  year-2  and  h  is  in  arcsec. 

The  distance  R  is  obviously  closely  related  to  the  semi¬ 
major  axis  of  the  orbit,  a.  If  we  write  R  =  a(R/a )  (for  rea¬ 
sons  that  will  become  apparent  shortly)  our  condition  for  de¬ 
tection  of  weak  curvatures  becomes 


47T2  (sin#)  M  T2  p 
24  [a  (?)]2 


(B2) 


The  strong  curvature  limit  is  defined  by  P  >  2 T,  where  P  is 
the  orbital  period,  given  by  P  =  a3/2  /  \[M,  where  M  is  the 
total  mass  of  the  system.  For  this  limit  we  therefore  have  the 
condition  at’!2  /\[M  >  2 T .  Rearranging  the  expressions  for 
the  weak  and  strong  curvature  limits,  we  obtain  the  limits  on 
the  semimajor  axis  that  define  the  fraction  of  stars  of  interest: 


(2 T)2/3  AI1/3  <  a  < 


Note  that  M  on  the  left  and  right  sides  can  have  different 
meanings:  on  the  left  it  is  always  the  total  mass  of  the  sys¬ 
tem;  on  the  right  it  is  the  total  mass  of  the  system  only  if  the 
motion  of  the  star  is  measured  with  respect  to  its  companion. 
If  the  star’s  motion  is  measured  with  respect  to  an  inertial 
frame,  then  M  on  the  right  is  the  mass  of  the  companion.  The 
expressions  therefore  can  be  used  for  investigating  planet  de¬ 
tection  by  the  astrometric  method,  in  the  long-period  limit, 
simply  by  setting  M  on  the  right  side  to  be  the  planet’s  pre¬ 
sumed  mass. 

We  have  chosen  e  =  2(e) /\/N,  which  is  twice  the  la 
uncertainty  of  any  angular  variable  derived  from  the  obser¬ 
vations  (assuming  N  ps  v).  Therefore,  the  quantity  p/e  ap¬ 
pearing  on  the  right  side  of  eq.  B3  can  be  thought  of  as  half 
the  signal-to-noise  ratio  of  the  parallax:  p/e  =  \p/ap.  (This 
holds  only  for  the  case  where  the  parallax  of  the  star  is  deter¬ 
mined  using  the  same  set  of  observations  used  for  the  proper 


motion  and  acceleration  determination.)  If  we  fix  the  p/ ap 
ratio,  eq.  B3  provides  the  semimajor  axis  limits  for  stars 
in  a  certain  distance  shell;  for  example,  if  we  wish  to  con¬ 
sider  stars  with  parallaxes  good  to  10%  ( p/ap  =  10)  when 
op  ps  10-3  arcsec,  we  will  obtain  the  limits  on  a  for  stars 
at  a  distance  of  approximately  100  pc.  It  is  interesting  to 
note  in  this  regard  that  if  we  have  a  parallax-limited  sam¬ 
ple  of  stars  uniformly  distributed  through  space,  the  average 
parallax  is  only  1 .5  times  the  minimum  parallax.  In  such  a 
case,  and  assuming  that  the  minimum  parallax  is  ap,  then 
p/e  <  1  for  most  stars  in  the  sample,  and  the  upper  limit 
on  a  will  generally  be  quite  restricted.  However,  most  cata¬ 
logs  are  magnitude-limited  and  their  stellar  distributions  are 
much  more  concentrated  toward  the  Sun.  In  such  catalogs, 
and  in  special  samples  of  nearby  stars,  there  will  be  enough 
stars  with  p/e  >  1  to  provide  a  useful  range  of  detectable 
semimajor  axes. 

To  assess  quantitatively  how  many  stars  fall  between  the 
limits  defined  by  eq.  B3,  we  need  to  know  something  about 
the  distributions  of  sin  6,  a,  M,  and  R/a.  Two  of  the  four  dis¬ 
tributions,  for  sin#  and  R/a,  are  easy  to  deal  with.  Although 
the  projection  factor  \/sin#  is  unknown  for  any  specific  bi¬ 
nary  system,  we  can  assume  that  the  direction  of  the  vector 
separating  the  two  components  is  randomly  distributed  over 
47 r  steradians.  For  such  a  distribution,  the  average  projection 
factor  is  0.87,  the  median  projection  is  0.71.  Therefore  we 
can  adopt  0.87  as  the  average  value  of  the  projection  factor 
%/ sin  #,  with  some  confidence  that  the  average  is  also  typical. 

The  distribution  of  R/a  values  presents  a  similar  situa¬ 
tion.  For  circular  orbits,  R/a=  1  at  all  times.  For  stars  in  ec¬ 
centric  orbits,  R/a varies  between  1-e  and  1+e,  whereeisthe 
orbital  eccentricity.  However,  binary  stars  spend  more  time 
near  apastron  than  periastron,  and  over  an  orbital  period  the 
average  value  of  R/a  is  1  t  /e1 .  Since  this  ratio  can  only  take 
on  values  between  1  and  1.5,  we  can  adopt  an  average  value 
of  R/a  of  1.3  without  much  concern  about  the  distribution  of 
eccentricities. 

Inserting  the  appropriate  numerical  quantities,  eq.  B3  be¬ 
comes 

1.59  T2/3  M1/3  <  a  <  0.56 /' (l/ -)  (B4) 

Since  e  =  2ap,  the  right  side  (upper  limit  to  a)  could  also 
be  expressed  as  0.61T(M p/op)1^2.  As  we  expect,  the  range 
of  applicable  a  increases  with  T.  For  M  =  1  and  p/e  =  1, 
7.4  <  a  <  8.7  if  T  =  10  and  34  <  a  <  87  for  T  =  100. 
These  limits  define  the  null  set  for  T  <  6.1  years  for  solar- 
mass  binaries  where  the  parallax  is  at  the  limit  of  detection, 
but  do  not  preclude  more  massive  or  closer  systems. 


